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ABSTRACT 

The  equation 

ut  =  (um)«,  i£R,  t  >  0 
u(-,0)  =  u0 

with  m  >  1  models  the  expansion  of  a  gas  or  liquid  with  initial  density  «o  in  a  one  dimensional 
porous  medium.  Denote  by  t  — ►  s±  (<)  the  vertical  boundaries  of  the  support  of  tt.  Caffarelli 
and  Friedman  have  shown  that  s±  €  C1(f±,oo)  where  t±  :=  sup{<  :  s±(t)  =  «±(0)}  is  the 
waiting  time.  Using  their  result  we  prove  that 


e±  e  C°°(t±,oo). 

Moreover,  we  show  that  the  pressure  v  :=  um_1  is  infinitely  differentiable  up  to  the  free 
boundaries  after  the  waiting  time.  Our  proof  is  based  on  a  priori  estimates  in  weighted 
norms  which  reflect  the  regularizing  effect  near  the  free  boundaries. 
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SIGNIFICANCE  AND  EXPLANATION 

The  equation  stated  in  the  abstract  describes  the  expansion  of  a  gas  in  a  one  dimen¬ 
sional  porous  medium.  While  existence  of  weak  solutions  can  be  obtained  by  standard  energy 
methods,  not  much  was  known  about  regularity  of  the  solution  near  the  free  boundaries.  The 
difficulty  is  that  the  equation  degenerates  at  the  boundary  and  the’^hyperbolic*  term  becomes 

dominant. 

ry-.-i  a  “hi  r 

We  show  in  this  report  that  despite  of  the  degeneracy  the  free  boundary  is  smooth  with 
the  possible  exception  of  a  discontinuity  in  the  derivative  at  th<?uSvaiting  time**.  Qtrr  method 
is  based  on  ^energy  estimates  in  weighted  norms  using  some  of  the  ideas  in  |6}.  -  V  -  *  ~  . 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies  with 
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C°°-REGULAR]TY  FOR  THE  POROUS  MEDIUM  EQUATION 
K.  Hollig  and  H.-O.  Kreiss 

1.  Introduction.  We  consider  the  porous  medium  equation 

ut  -  («m)«  =  o,  x  e  R,  t  >  0,  ^ 

u(-,0)  =  u0 

for  m  >  1  and  continuous  positive  initial  data  uo  with  connected  compact  support. 

It  is  well  known  [3,9, 10)  that  problem  (l)  has  a  unique  weak  solution  and  that  the  support 
of  u(-,  <)  remains  bounded  for  all  £,  i.e. 

suppu(-,f)  =  |r(t),  «(<)]. 

The  curves  r,  s  are  Lipschitz  continuous  |7],  but  in  general  not  Cl.  As  was  first  observed  by 
Aronson  [l]  r'  (and  similarly  s')  can  have  a  jump  for  t  equal  to 

tr  :=  sup{t :  r(t)  =  r(0)}. 

Caffarelli  and  Friedman  [4]  proved  that  a  classical  solution  of  problem  (1)  exists  up  to  the  free 
boundaries  for  t  >  max(tr,ta).  By  considering  the  equation  for  v  :=  ttm-1  (cf.  (2.1)  below) 
they  showed  that 

(i)  vtlvx,vxx  are  continuous  on  the  set  fl,  :=  {(x,t)  :  r(t)  <  x  <  s(t),  t  >  tr} 

(ii)  reC^tr.oo) 

(iii)  r'(t)  =  -^yu*(r(t),t),  <  >  U- 

The  corresponding  statement  holds  for  the  right  free  boundary  s.  In  particular,  the  functions 
in  (i)  are  continuous  on  the  closed  support  of  u  if 

®i(r(0))  V0(s(0))  *  0  (2) 

where  vq  :=  w(-,0).  With  the  aid  of  an  interesting  idea  of  Gurtin,  McCamy  and  Socolovsky  [5] 
it  has  been  recently  shown  [6]  that  r  E  Coo(0,7’j  if  vq  is  sufficiently  smooth,  (2)  holds  and  T 
is  sufficiently  small.  However,  this  method  does  not  yield  regularity  of  v. 

In  this  paper  we  obtain  the  following  optimal  regularity  result. 

Theorem,  v  €  C00^,.),  r  e  C°°(tr,  oo). 

Our  approach  is  different  from  the  method  in  (6);  it  is  based  on  the  smoothing  effect  of  the 
porous  medium  equation  in  a  neighborhood  of  the  free  boundaries.  We  prove  in  section  2  the 
following  a  priori  estimate. 
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Proposition  1.  Let  u  be  a  solution  of  (1)  for  which  v  G  C°°(fir)  and  assume  that 


«(0)  -  r(0)  <  k  1 

k  <  VqWO)),  |«ol  <  K~1  (3) 

l«o(r(°)  +  y)  -  »J(r(0))|  <  A(y),  y  <  k, 

where  k  is  a  positive  constant  and  A  is  a  smooth  function  with  A(0)  =  0,  A'  >  0.  Then,  for 
any  fc  G  N,  there  exist  positive  constants  6,  T,  A  such  that 

lrU,|T/2,T]  +  Mfc,n(tf,r)  <  A  (4) 

where  0 (6,T)  :=  {(*,<)  :  r(t)  <  x  <  r(t)  +  6,T/2  <  t  <  T }  and  |  |jt,n  denotes  the 
norm  on  W^( tl).  The  constants  S,T,  A  depend  on  k,  A,  fc;  in  addition,  T\A  depend  on 

|«>o|2fc+4,|r(0)+«/2,r(0)+<c]- 

In  section  3  we  show  existence  of  smooth  solutions  for  smooth  data. 

Proposition  2.  If  v0  £  C°°(suppv0)  and  (2)  holds,  then  v  G  C0O(suppt;)  and  r  G  C°°(0,  cx>). 

The  Theorem  follows  from  Propositions  1,2  by  an  approximation  argument.  Assume  that 
u  is  a  solution  of  problem  (l).  By  the  result  of  Caffarelli  and  Friedman,  (i)-(iii)  are  valid  for 
v  and  f.  Let  tf  <  t\  <  For  any  r  G  [<i,t2],  *>o  :=  v(-,r)  satisfies  the  assumptions  (3) 
of  Proposition  1  with  a  constant  k  and  a  modulus  of  continuity  A  which  depend  on  v,ti,t2 
but  not  on  r.  For  each  (fixed)  r  we  approximate  vq  by  a  sequence  of  smooth  functions  vo,j  £ 
<7°°(suppt/o)  for  which  (3)  remains  uniformly  valid  and  which  converge  to  v0  in  L00(suppuo). 
In  addition  we  require  that  (2)  holds  for  vo,j  and 

suppvo,,  =  suppvo 

v0ij(x)  >  0,  r(0)  <  x  <  s(0),  (5) 

suPluO,jl2*-M,[r(0)  +  */2,r(0)+«]  <  «>• 

J 

Let  denote  the  solutions  of  (1)  with  initial  data  uo  =  (vo,j)1^m_1^  By  Propo¬ 

sition  2,  Vj  G  C00(supptiJ).  Moreover,  the  conclusion  (4)  of  Proposition  1  is  valid  for  v3  and 
the  corresponding  left  free  boundary  Tj,  uniformly  in  j.  Passing  to  the  limit  j  — >  oo  it  follows 
that 

re^[r+r/2,r  +  r] 

v  G  ({(*.<)  :  r(t)  <  x  <  r(t)  +  6,  r  +  T/2<t<T  +  T}). 

Since  fc  G  N,  r  G  were  arbitrary  and  in  the  interior  of  suppv  the  regularity  is  known, 

the  Theorem  follows. 


2.  A  priori  estimates.  Troughout  this  section  we  assume  that  u  is  a  solution  of  (1.1) 
for  which  v  satisfies  the  assumptions  of  Proposition  1.  Substituting  u  =  in  (l.l)  we 

obtain 

vt  -  mvvxx  -  nv‘  =  0 
v(-,0)  =  v0 


-2- 


where  n  :=  l/(m  -  l).  The  change  of  variables 


y  =  x-r{t),  t >(*,<)  =  «>(y,t) 

transforms  the  left  free  boundary  to  the  vertical  axis  (y  =  0}.  Since  by  (ii 

Vt  =  -r'(i)  =  ntt)„(0,  t) 

the  problem  for  w  is 

wt  -  mwwyy  -  nujy  +  ntnv(0,  -)wy  =  0 
u;(-,0)  =  wo  :=  vQ( •  +  r(0)). 

For  the  proof  of  Proposition  1  it  is  sufficient  to  show  that 

l<?vti;lo,(o1,5)x|77217')  <  A',j<2k. 

We  need  several  auxiliary  Lemmas. 

Lemma  1.  /0#  /( y)2  dy  <  c,  y2(6~2f(y)2  +  f'(y)2)  dy. 

Proof.  By  scaling  we  may  assume  that  6  =  1.  Then, 

/V  =  /(!)’  -  2  /'»//' 

Jo  Jo 

<  /(l)1  +  1/2  J  I2  +  2  f 

where  the  first  term  on  the  right  hand  side  can  be  estimated  by  the  standard  Sobolev  inequality. 
Lemma  2.  sup0<,,<s  |y/(y)2|  <  e2  J*  y2{6~2f(y)2  +  /'(y)2)  dy. 

Proof.  Again,  by  scaling,  let  6  =  1.  Then, 

*/(*)’-/( l)2  -  f'  /(»)’  +  dy 

<  /o)1  +  if  r  +  j  v\n\ 

and  the  Lemma  follows  from  Lemma  1  and  the  standard  Sobolev  inequality. 

Lemma  3.  Let  Q(6,T)  :=  [0,£]  x  [0,T],  dQ  :=  [0, 6]  x  {0}  U  {6}  x  jO,T]  and  assume  that 
p  :=  minaQtw,,  >  0.  Then 

min  <  min  wv  <  max  wu  <  maxtv„. 
dQ  Q  Q  v  ~  dQ  v 

Proof.  Set  ri(t)  :=  ( p  -  t)exp(-it)  with  0  <  t  <  p.  We  differentiate  (2)  with  respect  to  y 
and  subtract  q'  4-  tq  =  0.  This  yields 


(2) 

(3) 


(u/„t  -  Vt]  +  \-mwwyyy j  +■  [((-m  -  2n)u;v  +  nu;v(0,‘))wvyJ  +  j  — c*y]  -  0. 


Assume  that  tvv(y,f)  =  r/(t)  where 

i  :=  sup{< :  >  rf(t)}. 

If  (y,  t)  €  Q\dQ  alJ  terms  in  square  brackets  are  nonpositive.  Since  y  /  0  this  is  not  possible, 
i.e.  we  must  have  q  <  tv„  on  Q.  Letting  <  -*0  proves  the  first  inequality  of  the  Lemma  and 
the  last  inequality  is  proved  similarly. 

Lemma  4.  If  26  <  k,  A(2£)  <  k/4 ,  then  there  exist  constants  T  and  C3  which  depend  on 
k  lvo|2fc-t--4,|^/2,*.]  such  that 

max  wv  -  min  wv  <  4A(6) 

Q(S,T)  Q((,T ) 

k/2  <  wv[y,t)  <  2 k'1,  (y,t)  €  Q(6,T), 

| d'/d$w{6,  t )  |  <  c3,  2i/  +  n  <  2k  +  3,  t  <  T. 


Proof.  The  maximum  principle  is  valid  for  problem  (1.1),  i.e.  u0  <  Uq  implies  that  u  <  u+ 
and  r~  >  r+.  By  (1.3)  and  our  assumption  on  6, 

v'0(y)  >  3«/4,  y-  r(0)  <  26. 

Using  this  and  (1.3), 

Vo  :=  max{0,  (y  -  r(0))(r(O)  +  26  -  y)/2}  <  v0  < 
max{0,  (y  -  r(0))(r(0)  +  4k~*  -  y)}  =:  v£ . 

For  the  solutions  of  (1.1)  with  initial  data  u*  =  (v*)1^m-1)  the  assertions  (i)-(iii)  are  valid 
with  tr  =  0.  Therefore,  by  the  above  comparison  principle, 

c  <  v(y,  t)  <  c_1 
-  c_,t  <  r(t)  -  r(0)  <  -ct 

if  6/2  <  y  <  36/2,  t  <  1.  The  constant  c  depends  on  6,  k.  We  choose  T'  <  1  so  that 

|r(t)  -  r (0) |  <  6/4 ,  t  <  T', 


which  also  yields 

c  <  t v(y,<)  <  c'1  if  36/4  <  y  <  5 6/4,  1  <  T'. 

On  the  rectangle  |36/4, 56/4]  x  [0,  T'\  the  problem  (2)  is  nondegenerate  and  the  last  inequality 
in  (4)  follows  from  parabolic  regularity  theory  if  T  <  T'  [8] .  We  set  T  :=  minlT7,  A(£)/c3}. 
Then 

K(M)  -  <  ^K,(M")i  <  m 

C3 

which  yields  the  first  two  inequalities  for  (y,  <)  €  dQ  and  therefore,  in  view  of  Lemma  3,  also 
for  (y,  <)  €  Q. 


% 


I 


f 
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Proof  of  Proposition  1.  Let  0  =  T_i  <  To  <  ...  <  Tu+i  =  T/ 2.  We  prove  by  induction 
on  l  that  for  sufficiently  small  6, 

^  T  b 

max  f  ydi,+ 1  u;(y, f)2  dy  +  f  [  y2d‘‘t'2tn(y,t)2  dydt  <  A"(f),  0</<2fc+l.  (5) 

Ti<t<T  j0  jTi  j0  " 

The  constants  A"  depend  on  k,6,  X,k,Tv,  |uo!2*+4,|/>/2,ic)-  By  Lemma  1, 

|ajti>(-4t)*|o,|o/i  <  eg  J  dlw(-,t)2  +  dJv+'w(-,t)2 

<  cscx6~2{A"(j  -  1)  +  2A"{])  +  A"(]  +  1)) 

which  shows  that  (5)  implies  (3). 

Since  w  and  wy  are  bounded,  inequality  (5)  is  obviously  valid  for  /  =  -1.  We  assume  that 
(5)  holds  for  l  <  j  and  set  W/(y,t)  :=  d3y+1w(y ,f  +  T,_i).  Differentiating  (2)  (j  +  1)  times 
with  respect  to  y  and  replacing  t  by  t  +  Tj-\  we  obtain 

(W,)t  -  mW-\Wj+2  ~  ((2»+  (i  +  l)m)W0  -  nW0(0,  •))%  +  ,  -  £  cvllWvW»  =  0  (6) 

V+*=J+ 1 


where  are  constants  which  depend  on  j.  We  multiply  (6)  by  t2yWj  and  integrate  over  the 
interval  [0,  $], 

\(ft2yW*  dy)t  +  m  J*  t'yW-iW?^  d y  = 

/ tyW > 

+  mt36W-l(6,t)Wj+i(6,t)Wj{6,t) 

r  (7) 

-m  / 

+  J  t2y\(2n  +  (j  +  l)m)W0  - 

+  f  t2yWvW„Wj. 

The  third  term  on  the  right  hand  side  of  (7)  equals 

-mt2(W_,(M)  +  Wo(M))W;(M)2/2  +  m  J  t2{W0  +  yWx/2)W2. 

Proceeding  similarly  with  the  fourth  term  on  the  right  hand  side  and  using  (1.3)  and  (4)  we 
deduce  from  (7)  that 


i  </>*■/>. +  t7/v,v-  s 

C4C3  +  J  tyW2 

-  J t2\-mW0.+  (n  +  (j  +  l)m/2)W0  -  ^W0(0,  -)}W2 

+  ««  ,<*!£?<,  I  / 

«'+#  =  >+!  ^ 


-5- 


where  the  constant  c*  depends  on  k  and  the  constant  eg  depends  on  j.  We  estimate  each  of 
the  integrals  appearing  on  the  right  hand  side  of  (8)  seperately.  By  the  definition  of  Wi  and 
the  induction  hypothesis 


/  tyW]{y^)2  dy  < 

J°  (9) 
«  J  t2Wj(y,t)2d  y  +  J  y2d3y^lw(y,t  +  Tj-i)2  dy. 

By  (4),  |W0(y,  t)  -  W0(0,t)j  <  4A(6)  and  /c/2  <  W0(0,t)  <  2/c-1.  Therefore  the  term  in  square 
brackets  in  the  second  integral  on  the  right  hand  side  of  (8)  can  be  estimated  by 

—  Cg,  if  j  =  0 

n/c/4  -  CgA(6),  if  j  >  0  (io) 

>  n/c/4  -  C6 A(^ )  -  max(0, 1  -  j)c6 

where  c6  depends  on  j,/c.  Finally  we  estimate  |  /  ^yW^W^Wj],  Set  W0(y,t)  :=  W0(y,t)  - 
Wo(0,t).  Integrating  by  parts  and  using  (4)  it  follows  that 

I J  <2y*Vi%2|  < 

t26W0{6,t)Wj{6,t)2\  +  \  J  t2W0W2\  +  2|  J  t2yW0WjWj+\\  <  (11) 

4\(S)c2  +  8\(6)  J  t2W2  +  8\{6)  J  t2y2W2+l 

if  6,  t  <  1.  We  have 

|  J  SyW^W^W, |  <  c  j  +  (12) 

vhere  Bytl(t)  :=  f  t2y2W2W2.  If  v  <  fi  <  j  it  follows  from  Lemma  2  that 

£  <2(rjmaxf  x  (j  yW„(y,t)!  4  y) 

<  «»/-*( J y2(w?  +  wy\,))  x  ( J vK)- 

Therefore,  using  the  induction  hypothesis, 

B^{t)  dt  <  c26~2(A"(v  -  1)  +  >1»)  x  A"(»)  <  c7A "(j  -  l)2.  (13) 


-6- 


Combining  the  estimates  (9-12)  it  follows  from  (8)  that 


|  (/*»?),  +  ™  /  tV«i?+1  < 
c4e|  +  £  J  t2W2  +  i~lb{t) 

-  {m c/4  -  c6A(6)  -  max(0, 1  -  j')c6)  J  t2W 2 
+  c5(4A(6)c2  +  8A(6)  J  t2W2  +  8A(«)  f  t2y2W2^) 
+  c6(£  j  t2W2  +  r'  4  max  . 


where  6(t)  =  /  y2dj+1u>(y,t  +  T^-j)2  dy.  We  choose  6,«  so  that 


.  «  TTO/C 

8c5A(6)  <  — 

4 

£  +  CeA(6)  +  8cbA(^)  4-  c$£  <  tik/4. 


Then  we  obtain  from  (14)  that 


i(/‘V*7),  +  £ 

c4c|  +  £*J6(<) 

+  c6  max(0, 1  -  ji)  /  ^ 

+  C5£_1  maxBv<1(t). 


Since,  induction  hypothesis, 


T—T 

f  b(t)  dt  <  A"(>  -  1) 

Jo 


it  follows  from  (4)  and  (13)  that  for  any  t  €  [0,  T  -  Ty_ij, 

If  t*yWAy,t)2  dy  +  x/  /  Ty2^>+i(y^)2dydr  < 
c4c^  +  £_1A"(j  -  1)  +  4t3K~2  +  chr\1A"(j  -  l)2t. 


This  completes  the  induction  step. 


3.  Existence  of  smooth  solutions.  In  this  section  we  outline  the  proof  of  Proposition 
lich  justifies  the  approximation  argument  in  the  introduction.  Similarly  as  in  section  2 
ransform  the  equation  (2.1)  to  a  fixed  domain.  Let  £  6  C°°[0,  lj  satisfy  £'  <  0,  0  <  £  < 
(y)  =  1  for  0  <  y  <  k,  £(y)  =  0  for  2*c  <  y  <  1  and  set  »j(y)  :=  £(l  -  y).  Assuming 
out  loss  that  r(0)  =  0,  s(0)  =  1  the  change  of  variables 

y  =  *  -  £(y)r(0  -  v(y){s{t)  -  i) 

t >(x,t)  =  w(y,t) 

sforms  the  free  boundaries  to  the  vertical  lines  {y  =  0}  and  {y  =  1}.  One  easily  verifies 
the  transformed  equation  for  w  is 

“’t  -  {m/x2)wwyy  -  (n/ x2)wl  +  (n/x)£wv(0,-)wy  4-  (n/x)»?twv(l, -)wv 

+  (mXy/x3)«>u’y  =  0,  0  <  y  <  1,  t  >  0,  (2) 

u>(-.0)  =  w o  :=  t;0 


e  ^ 

x(y^)  =  1  -  n£'(y)  f  tWy(0,r)  dr  -  nr)'(y)  j  u>„(1,t)  dr. 

Jo  Jo 

neighborhood  of  the  left  boundary  {y  =  0}  we  have  x(y)  =  1  and  equation  (2)  coincides 
equation  (2.2).  Therefore  an  analogous  a  priori  estimate  is  valid. 

ma  5.  Assume  that  w  6  C°°([0, 1]  x  [0,  T])  and  that  Wq(0)  Wq(1)  /  0.  Then  for  any 
l 


(  max 

vo  <t<T 


<xr/  y(l~y)d iw(y«Oady)  +  (fQ  fo  y2(i  -  y)2^t+,ty(y,02  dydt)  < 


(3) 


e  c  depends  on  l,T,v q. 

proof  of  this  Lemma  is  completely  analogous  to  the  proof  of  Proposition  1.  Instead 
ultiplying  equation  (2.6)  by  t2yWj,  we  multiply  the  corresponding  equation  obtained  by 
entiating  (2)  by  y(l  -  y)d3+lw(y,  t).  Because  of  the  weight  y(l  -  y)  no  boundary  terms 
ar  when  the  appropriate  terms  are  integrated  by  parts.  The  estimates  are  somewhat  more 
dicated  because  of  additional  terms  involving  x •  But,  these  complications  are  merely  of 
lical  nature. 

Given  the  above  a  priori  estimate  it  is  straightforward  to  prove  a  corresponding  local  ex- 
ce  result  via  finite  difference  or  finite  element  approximation.  This  completes  the  (outline 
e)  proof  of  Proposition  2. 


f 


'  r 
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